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We have developed a new approach based on matrix product representations ol ground states 
to study Quantum Phase Transitions (QPT). As confirmation of the power of our approach we 
have analytically analyzed the XXZ spin-one chain with uniaxial single-ion-type anisotropy and our 
results compare favourably with previous numerical studies. In addition, our description lets to 
know which part of parameters space of the Hamiltonian is most likely to be exactly solvable. 
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PACS numbers: 

The challenging problem of determining the ground 
state of a quantum spin system and its properties such 
as the energy and degeneracy, and the behavior of asso- 
ciated correlation functions and QPTs [H developed sev- 
eral numerical and analytical approaches in many-body 
physics. Numerically, the density matrix rcnormalization 
group (DMRG) @f and analytically finitely correlated 
states [3|, also known as translational invariant matrix 
product (TIMP) states 0, 110,0,11 as generalization of 
valence bond states @ and several other methods pro- 
vide a valuable framework to deal with quantum spin 
chains, although still there is lack of specific formalism 
for studying QPT points. 

Recently, in [Io| a method was proposed to study a new 
family of QPTs based on the MPS formalism. They de- 
rived the correlation functions in terms of the eigenvalues 
of the transfer matrix E, and have shown that the cor- 
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diverges when level crossing 



occurs between two largest eigenvalues of E, v\ and v^. 
This is an interesting kind of QPTs which can be found 
for MP states [H EH, EH- However, it is clear that in 
such QPTs there is no change in the ground state energy 
from one side of the critical point into the other side, 
because the ground state energy is always zero, which is 
far from a standard expectation for physically relevant 
QPTs. 

In this letter, we propose a method to characterize 
QPTs based on MP representation of the ground state 
and we support our general discussions by comparing 
our analytic results with the numerical results of a well- 
known model for spin one chains. We also propose a for- 
malism to find an MP representation for ground state of 
a given Hamiltonian, and present a criterion to determine 
in which part of the parameter space of the Hamiltonian 
the ground state is not an MPS. Based on this descrip- 
tion we show that subspaces accompanied by the higher 
degeneracy in the local ground state are the most likely 
to yield the ground state exactly. This gives a new idea 
to find exact ground states of several models, and, using 
other methods like spin waves to determine ground state 
properties near the exactly solvable subspace. 

We can construct the parent Hamiltonian for |^) = 
^J2t.. iN ^(Ai 1 ...Ai N )\h...i N ), by finding fc« = 
Sj =1 Oi|ui)(ui|, which introduces the interaction among 



k particles in a translational invariant parent Hamilto- 
nian H = Y,f =l hi._i + ki where the \viYs span the null 
space of the reduced density matrix of k neighbor- 
ing sites. For at > 0, H is positive semidefinite and |\&) 
is the ground state with zero eigenvalue. Constructing 
the parent Hamiltonian for an MPS is based on the fact 
that the interaction among the particles can be modeled 
in a finite-dimensional vector space (which comes from 
the definition of the finitely correlated states). Despite 
this, the constructed parent Hamiltonian is typically not 
similar to well-known models, for example the Hciscn- 
berg model. So, it is important to study a formalism 
to find MP representations for the ground states of arbi- 
trary models, and to aim at understanding in which part 
of the parameters space of a given Hamiltonian an MP 
representation for the ground state is possible. 

We begin with the Hamiltonian of a translationally 
invariant local interaction, 
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It is possible to represent the local interaction in terms 
of its eigenvalues and eigenvectors by spectral decompo- 
sition, = J2i=i a i\ e i)i e i\, and also, to determine the 
smallest eigenvalue of h^ k \ Let us consider the given 
Hamiltonian, with periodic boundary condition, as the 
sum of a positive semidefinite Hamiltonian and a con- 
stant term: 
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where a is the smallest eigenvalue of the and H' is 
a positive semidefinite operator, because 



fc'(fc)= s £ r »( a ._ a )| e .)< e .| 



(3) 



where g is degeneracy of a and a* > a. 

Now, we pose the question: "Is there any set of matri- 
ces {Ai} to construct the MP |^) such that = 



2 



0?" That is, can one arrange things such that the re- 
duced density matrix of \*f?) for k consecutive sites has 
|ej)s as the basis for its null space? The MPS formalism 
indicates that the |ej)s impose a set of constraints on the 
AiS, so this enables us to look for d matrices to repre- 
sent MP ground state for H' with zero eigenvalue. As we 
now discuss, thinking along these lines yields a criterion 
- Theorem 13 of Q- to know whether there is a TIMPS 
as ground state. 

It is clear that there is a bound on the ground 
state energy E g of H for N particles, namely E g > 
Not. Assuming that there exists an MPS such that 
{^mps\H'\^mps) = 0, it can be shown that \^mps) 
is eigenstate of H with E g = Na. Hence, if there 
exists {Ai}f =1 that satisfy the equations imposed by 

the {|ei)} i=1 , then the ground state energy will be 
Eg = Na. The greater the degeneracy in the local 
ground state impose the smaller the number (d k — g) 
of constraints on the AiS. Intuitively, the more degen- 
erate subspaces in the space of parameters {</} are ones 
for which it is more probable to find an exact ground 
state [HQ- For example, AKLT have solved the H = 
E"i(<Si-$+i) -P(Si.S i+1 ) 2 for = -|, exactly. It can 
easily be shown that (3 = — ^ is the crossing point of the 
non-degenerate local ground state of energy —2 — 4/3 (for 
(3 > — |) with the triply degenerate local ground state of 
energy —1 — /3 (for (3 < — |). So, at the level crossing 
point (j9 = — 5) the local ground state degeneracy is max- 
imum (g = 4) , and it can be shown that the local Hamil- 
tonian consists of the 3 2 — 4 = 5 remaining vectors which 
equal the spin-2 projector P- i+1 = J2 m =-2 I 2,m)(2, m |. 
Fortunately, because of the continuity of E g in terms of 
the parameters of the Hamiltonian, i.e. coupling con- 
stants {J}, Eg — Na approaches zero, smoothly: 

lim E„-Na = 0, (4) 
{J}^{J}„ 

where, {J}b identifies the boundary of the region in which 
the ground state has TIMP representation. This inves- 
tigation can be also performed using the DMRG to 
calculate E g . These facts allows us to look for regions 
in the parameter space of H, where the ground state is 
a TIMPS. When the ground state energy E g is equal to 
the minimum possible energy Na, the ground state mini- 
mizes energy locally. This means thatTr(/i( fc V fc )) = a or 
Tr^'WpW) = 0, where ftw is the local Hamiltonian for 
k neighboring sites, h'^ is the positive local Hamiltonian 
in Q and p^ k > is the reduced density matrix for k sites. 
In another language tr (h'^'p^n = means that corre- 
lations between k particles can be specified by a finite 
dimensional vector space. The equality E g = Na leads 
to tr(h'^ p^) = which means that h'V°> can be written 
in terms of the null space of p^ 3, i.e., the eigenvec- 
tors of the local Hamiltonian specify the null space of the 
reduced density matrix for each fc-site block. So we can 
expect that, in these cases, the global ground state can 
be determined by studying local ground states like the 
MPS. 



Taking a step further on ([3]) for a given H, we now 
introduce a method to study QPTs. Consider the de- 
composition of hi . It is clear that the eigenvectors 
corresponding to the smallest eigenvalue were omitted. 
Each remaining eigenvector results in a constraint equa- 
tion on the matrices AiS. For example, a vector like 
|e) = (l,a,b,c), which comes from a two body Hamil- 
tonian for spin one- half particles, imposes an equation 
like A + A + + aA + A_ + bA_A + + cA_A_ = 0, where 
A± denote matrices related to |±) states of spin one-half 
objects. 

Therefore, a level crossing in the smallest eigenvalues 
of h leads to at least one change in the eigenvectors that 
describe h' , and so a change in the set of equations gov- 
erning the AiS. It is possible to have more alternations 
when a is degenerate. Solving these equations for a spe- 
cific D, which is the dimension of A^s, yields different MP 
ground states for each individual region, because they 
were obtained from different sets of equations. As such 
we can potentially identify a QPT purely by examining 
the structure of the low energy states of the h. 

This is only a rough idea to give insight into our pro- 
posal. It is important to relate these rough ideas to dis- 
continuity in a physical observable or define an order pa- 
rameter to determine QPTs by this description. 

Consider a model which has at least one level crossing 
in its smallest eigenvalue in terms of the Hamiltonian pa- 
rameters. Thus, there are two corresponding regions and 
in each region, one of eigenvalues of h is the smallest. At 
the level crossing point (it can be a point, a line, etc.) 
both of them are the smallest eigenvalues. Without try- 
ing to solve the equations in both regions, using Eq. ([¥]), it 
is possible to find out whether the ground state is Trans- 
lational Invariant MPS or not. The cases in which the 
ground state in different regions arc TIMPS, looking for 
discontinuity in derivatives of the ground state energy 
enables us to study QPTs. For instance, if a and j3 are 
the smallest eigenvalues of h and if the numerical inves- 
tigations show that Eg = Na and E g = Nf3 in the two 
regions, any discontinuity in derivatives of E g - or non- 
zero derivatives of a — (3 in terms of {J} - implies that a 
QPT occurs in the cross point. 

There is a possibility of finding an order parameter 
in order to characterize QPTs in another way. Let us 
start with the example of a spin one-half Hamiltonian 
with a two body interaction. Consider a model that pos- 
sesses \e\) = (1,0,0,0) as eigenvector corresponding to 
the smallest eigenvalue a in one region, but are higher 
energy vector for the other region. This means that in the 
first region A + A + ^ and in the other region A + A + = 0, 
because a is not the smallest eigenvalue there. So, due 
to I , in the second region, the coefficient of all terms 
with consecutive up (|+)) states are zero. In this case 
O =| ++)(++ I is a suitable local operator to define 
an order parameter, because it can be shown that the 
quantity lirn/v— >oo S^ =1 < 0,.i+i > is non-zero in the first 
region and zero in the other one. 

Lets keep considering two body Hamiltonians, which in 
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practise are the most physically interesting. The equa- 
tions resulting from the spectral decomposition of the lo- 
cal Hamiltonian are linear combinations of the products 
AiAj. Clearly, these equations are nonlinear in terms of 
the j4,'s, and the parameters in them. But, it is possible 
to easily linearize these equations by replacing AiAj with 
matrix My. Again, consider that a level crossing occurs 
for the smallest eigenvalue of a local Hamiltonian. Now 
there are two sets of equations corresponding to the two 
possible regions in the space of parameters. The equa- 
tions governing the MyS can be solved, which leads to 
the relation between My-s in both regions (after lineariz- 
ing the equations). There are always solutions to the 
Mij's, however this does not mean that there are AiS 
which satisfy all equations, because in the regions where 
Eg > Na solutions are accompanied by a contradiction 
between My = AiAj and My-s. 

Now, based on the above example and discussions, we 
see that it is possible to define an order parameter by 
studying solutions of linearized sets of equations in mod- 
els that possess a level crossing in their smallest eigen- 
value. In general, 
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corresponds to an operator O a p = \aj3){aj3\ which can 
be an order parameter when the sets of equations imply 
that 3a, (3 3 {1, ...,d} such that M a p = Mp a = in one 
region and M a p = B a Bp ^ in another region. Clearly 
the above order parameter is zero in the region where 
M a( 3^o and is non-zero in another region. In other words, 
there is no ordering like | af3a...) in one region while there 
is such a pattern in another region, even when there is 
no TIMP ground state in one of the regions. 

In summary, the above condition yields a method to 
detect a QPT. This condition, based on the formalism for 
finding MP representations of ground states, lets us study 
non-trivial QPTs in some models. In the next section 
we study a well-known model for spin one systems and 
compare our analytic results with numerical outcomes. 

We now apply our previous discussions on a well-known 
model, and analytically derive the presence of a QPT 
which previously was only known to exist through numer- 
ically studies. Consider the following two body Hamil- 
tonian as the parent Hamiltonian for a sp in one chain, 
which has been studied numerically in [15| . 



yields eigenvalues and eigenvectors. The eigenvalues are 

( ex = J Z + D 
e 2 = {-J, + D + y JJl - 2J Z D + £> 2 + 8 )/2 
e 3 = (- J z + D- y/ Jl - 2J Z D + D 2 + 8)/2. 



e 4 = -J z + D 



e 5 = (D + VWTl )/2 
[e 6 = (D- y/W+l)/2 



A simple calculation shows that there is a level crossing 
for the smallest eigenvalue, which is between e\ and e 3 . 
This means that there are two regions in the space of 
parameters {J Z ,D}. The following equation shows the 
line that separates these regions, which can be derived 
from ei — e 3 = 



J z = (-D - 

Let us denote regions R\ 
and R 2 : J z < {—D - 



V# 2 +4)/2 
> ( 
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In Ri the 



y/W+A)/2 

smallest eigenvalue is e% and the corresponding eigen- 
vector is (0, 0, 1, 0, — e 2 , 0, 1, 0, 0), whereas in R 2 the 
smallest eigenvalue is e\ with (1,0,0,0,0,0,0,0,0) and 
(0, 0, 0, 0, 0, 0, 0, 0, 1) as eigenvectors. 

Now, it is possible to look for TIMP ground state in 
each region by solving the corresponding sets of equations 
for each region. Omitting the above eigenvectors and 
translating the remained ones leads to the following sets 
of equations: 



i?! : 



Mn = 
Mil = 

Mn - e 3 M 00 + M Tl = 
Mn - Mn = 



e 6 M i 
e 5 Mi 
e 5 M i 



M I0 
M i 
M 
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e 6 Mi + M 



Mn " e 3 M 00 + M n = 
Mn - e 2 M 00 + M lx = 
Mu - Mn = 
R 2 : e 6 M i - M l0 = 
e 5 Mi + Moi = 
e 5 M i - M I0 = 
e 6 Mi + Mqi = 
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N FIG. 1: E a — (N — l)a in terms J z for an open chain for 

H = V S l x S l x +1 + S l y Sl +1 + J z S l z Sl +1 + D(Sl) 2 (6) N = 30 and D = 0. For J z > -1, in Ri, a = e 3 and for 
,■_ -i Jz < — 1, in i?2, a — ei. 



A spectral decomposition of the local Hamiltonian easily 



Numerical results for N = 30 particles, shown in Fig. 
Q] indicate that in R 2 there is a TIMP representation 
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FIG. 2: The solid line is Eq. © and is the result of our 
consideration based on changes in the structure of equations 
governing matrix product representations of the ground state. 
The other parts of the figure come from numerical studies 
performed in JT^]. 

for the ground state, whereas in R\ DMRG shows that 
condition Eg = (N — l)a is violated. This means that 
the ground state is not of the TIMP form, so there is 
no solution for the equations governing the matrices to 
represent a TIMP ground state. Note that the numerical 
simulations using DMRG were performed for an open 
chain which clearly consists of an N — 1 block so the N 
in ([4} must be replaced by ( N — 1) . 

Since E g — (N — l)a does not vanish in the ground 
state is a TIMPS only in i?2, and there is no solution with 
finite dimensional AiS in R\. Solving the set of equations 
governing the AiS in i?2 shows that only Mu = A\ A\ and 
Mil = A\Ai are non-zero. This means that i?2 is the fer- 
romagnetic phase of this model. It is clear that M\\ and 
Mu are zer0 m -^2, so local operators On = |11)(11| 
and On — 1 11) (11 1 " an< i an y linear combination of them 
- are suitable choices to show that a QPT occurs in the 
boundary ([7]) between R\ and i?2- In this case one can 
select the more familiar operator S\. Hence, the line ([7]) 



demonstrates a boundary between a ferromagnetic phase 
and other (generally anti-ferromagnetic) phases, of this 
model. Figure [2] shows our analytical result in compar- 
ison to numerical studies using exact diagonalization in 
[IH which is a good evidence for confidentiality of our 
results. 

We have proposed a method for studying spin chains 
with local interactions which in certain circumstances al- 
lows one to identify regions in parameter space for which 
there is no TIMP representation of the ground state. In 
addition, we have shown that in the space of parame- 
ters of the local Hamiltonian, subspaces where the lo- 
cal ground state has the highest degree of degeneracy 
are the most likely to be exactly solvable which can be 
guiding in further studies. Furthermore, it is possible to 
use this description for identifying and studying quantum 
phase transitions. In contrast to a family of QPTs that 
have been studied in [n|, QPTs that can be considered 
through this description are accompanied by a change in 
energy. 

We can use this formalism to study one dimensional 
spin systems and to benefit from the simplicity of the 
Matrix Product representation to calculate correlation 
functions, correlation length or investigate entanglement 
in these systems. In addition, the possibility of finding 
an approximate solution for the set of equations, i.e., a 
TIMPS which approximately describes the ground state 
for example by Ai& satisfying tr(J^ - =1 d CijAiAj) = e 
for a small finite e rather than e = can be studied. 

Finally we mention that it is possible to turn things 
around and impose level crossings of the smallest eigen- 
value to construct and engineer new models with a 
TIMP ground state accompanied by determined Quan- 
tum Phase Transition, and study behavior of correlations 
and entanglement near critical points. 
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